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Abstract. We model the effects of self-excitation/damping and shock transmission of Alfven 
waves on solar-energetic-particle (SEP) acceleration at a coronal-mass-ejection (CME) driven par- 
allel shock. SEP-excited outward upstream waves speedily bootstrap acceleration. Shock transmis- 
sion further raises the SEP-excited wave intensities at high wavenumbers but lowers them at low • 
wavenumbers through wavenumber shift. Downstream, SEP excitation of inward waves and damp- 
ing of outward waves tend to slow acceleration. Nevertheless, > 2000 km/s parallel shocks at ~ 3.5 
solar radii can accelerate SEPs to 100 MeV in < 5 minutes. 
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INTRODUCTION 

Self-excited Alfven wave turbulence near CME-driven shocks enhance SEP scattering 
and first-order Fermi acceleration [1, 2]. Steady-state models [e.g. 3, 4], which do 
not address acceleration time scale, consider only self-excitation of upstream (US) 
waves, because the constant-density SEPs downstream (DS) are purely convected. Time- 
dependent models, however, must treat DS wave evolution as well, because the SEPs 
stream down density gradients in both directions away from shock. For simplicity we 
consider energetic protons interacting resonantly with circularly polarized parallel and 
anti-parallel Alfven wave modes (R^, L^) near a parallel shock and focus on two 
processes affecting wave evolution: SEP excitation/damping and shock transmission. 

We find that US, protons accelerated up to energy E amplify outward waves and damp 
inward waves down to some minimum wavenumbers k. Each SEP- amplified US Alfven 
wave mode then crosses the shock to emerge DS as a small inward mode plus a dominant 
outward mode, amplified and shifted to higher k by shock transmission, consistent with 
[5, 6, 7, 8]. DS, first, the /^-shifted waves resonate with lower-energy protons. Secondly, 
SEPs scattered back to the shock by inward waves lose energy compared to those 
scattered by outward waves. Thirdly, the inward streaming SEPs damp outward waves 
and amplify inward waves. 1 These DS processes slow acceleration and soften the SEP 
energy spectrum at high energy, consistent with the finding in a simpler model by [10]. 
We present model calculations below to quantify these effects and the influence of the 
outer boundary condition. 


1 Inward/outward relative to the plasma. All modes travel inward everywhere relative to the shock. 



MODEL 


We adopt a radial mean magnetic field B = Bq (ro / r ) 2 and plasma proton number density 
n K — n H,o( r o/ r ) 2 ’ h^ n ce Alfven speed — UA,o( r o/ r )> with r the radial distance. 
.Bo = 0.143 gauss, ^h, 0 — 2 x 10 5 cm -3 , and Va, 0 = 700 km s _1 are typical values 
at ro = 3.5 r 0 [9, 11]. Relative to the Sun, the shock, US, and DS plasma travel at ^sh = 
2500 km s 1 , V u = 83 km s , and V d = 1880 km s , respectively. At time t — 0, the 
shock is at r S h = ro, fluid compression ratio Cf = 3.895, and Mach number M& = 3.45. 
The outward and inward Alfven wave velocities are W±(pt) = V(r.t) ± V A (pt). The 
SEP and Alfven wave distributions evolve according to 

d,F + d r [Qiv + W+)F] + (fiF) + Pdp(FP/P) = 

dfi(DfifidfiF +DfipdpF) +P~ 2 dp[P 2 (Dp^d^,F +DppdpF)\, (1) 

fi ={v + flW+)( 1 -F 2 )/r- [(1 -M 2 )/v] {d t W+ + (nv + W + )d r W + }, (2) 

P/P = -W + (l-ii 2 )/r- (ju/v) [d t W + + (nv + W+)d r W + }, (3) 

d t V a +d r (W a '¥ a )+r,d n {(2W a /r-dW ( y/dr)'¥ a } = y a ^ a , (4) 

D m = (y/2PfY J dkI a R%, (5) 

<7 

7a (*, r,t) = 2nW go V K f jdlidPP^R^Gf /[S 2 {\ - pV a /v) 2 }. (6) 

In the above, F(P,fi,pt ) = f (P. ii.pt) BoP / (BP / ) , with / the SEP phase-space den- 
sity in mixed coordinates: r, t in inertial frame; velocity v, rigidity P, and pitch cosine 
jU in local outward wave frame. Po = a constant. D a p = momentum diffusion tensor 
[12]. Besides focusing and adiabatic deceleration, fi and P in eqs. (2) and (3) include 
differential frame transformation of (/J, P) following particles. I a (k.pt) — magnetic in- 
tensity of a- mode Alfven wave. y ¥ a ( / q,pt) = 2I a r\\W a \/ (tJoVa) 00 wave action density, 
T] = k/B, rjo = a constant, and y a = fractional growth rate of I a . In eq. (6), S = total 
particle energy; Gf = d^f - (jU<^/ - Pd P f)V a /v; R° il (ji.P.k.V A ,B) = wave-particle 
resonance function [12]; g a = 1 , V a = 0 for a = P+,P+; and g a — —1, = —2V A for 

< 7 = R—,L 

US, I R+ = I L+ oc k -5 / 3 are initialized to a steady-state solution of the wave kinetic 
equation (4) and I R ~ = I L ~ = 0.04 I R+ \ these are normalized to give mean free path 
A > 0.4 AU at 1 MeV [12]. DS, P 7 and n\\ are derived via shock transition. At t = 0, the 
DS shock-transmitted waves are already enhanced by a factor of « 80. We specify US 
and DS isotropic seed protons / seed = 5bnu(r sh )P~ & /{4 k(P- 5 -P fe “ 5 )], (P a <P < Pb ), 
with b = 0.0028, P a = 5.92 MV (19 keV), and P b = 9.95 MV (53 keV). 

Eqs. (1) and (4) are recast in x where x = r — r S h(f) (hence V^, etc. become time- 
dependent) and solved numerically via finite difference approximation using operator 
splitting and dynamic time steps [12] in a box comoving with the shock. There are 200 
DS cells and 400 or 2000 US cells. The grid is defined as follows: xj — xj + yAx, j = 
-200:400 or 2000, Ax = 2. 14 x 10“ 4 r G ; /z, : = (/+ l/2)A/i., i = -20:19, A ji = 0.05; In P £ = 



In 1)A InP, 1=1 :46, A InP = 0.08664, P x = 10.38 MV; and log77 m = log t]i + 

(m — l)Alogrj)m= 1 : 48, A log 7] = 0.05, logr/i = 4.87 x 10 -4 MV -1 . 

At the moving inner and outer boundaries, we specify 201 = Fu , j =- 200. 

Fi,ej=m = Fi/j= 400. 'P(T,mj=40i = ^ambient- Protons are accelerated via P-transport 
in cells adjacent to the shock (eq. (3)), where \ixAW + /v£\ < 0.2AlnP; otherwise, frame 
transformation is performed on shock-crossing protons to evaluate more accurately the 
j u and P increments. Each US wave mode IJJc u ) crossing the shock is converted pre- 
serving helicity into two DS modes Id{kd) as follows: [5, 6], 

Id(k d )/Iu{ku ) = (l+Ud/VcJ) 2 c w /4, (7) 

C* = Cf(M A -Zu)/(M A -Ct/y/Cf ), (8) 

c w = kd/ku = wave compression ratio, £„ = ±1 and Q = ±1 for outward/inward US 
and DS waves. Since nonlinear processes limit wave growth, we impose l c < l sa t = 
B 2 / {'ink) so that X > 3 r g , r g = gyroradius, to stay within quasilinear theory. 


RESULTS 

Figures la and lb show the coupled evolution of the proton intensity Je versus E and the 
Alfven wave intensity I R+ versus k/B, respectively, just US of the shock in the reference 
case. In 5 minutes protons are accelerated to 100 MeV and wave growth spreads from 
k/B « 0.3 MV -1 to 0.002 MV -1 . 20 - 50 keV seed protons emerge across the shock 
anisotropically and excite ambient R+ and L+ Alfven waves at resonant wavenumbers 
k ~ B/[P(ji — V a /v)\. The amplified waves enhance scattering, providing a positive 
feedback loop that rapidly bootstraps acceleration and wave growth. The strong wave 
growth results in the wave spectra being limited by I sat . 

Snapshots of the radial profiles of jg and I R+ are shown in Figs, lc, Id for 1.3 MeV 
protons and le, If for 68 MeV protons, respectively. At 1.3 MeV, proton acceleration 
“ignites” at ~ 18 s and is in full throttle by 120 s. Simultaneously, the US X falls from 
0.4 AU to ~ 3 x 10 1 AU (Figs, lc and Id). In contrast, acceleration at 68 MeV ignites 
later at ~ 150 s, when X falls below 1 x 10 -4 AU (Fig. If), and does not reach full 
throttle at 300 s (Fig. le). 


TABLE 1 . Model Runs 



Upstream 

cells 

D^p.Dpp 
X= 0 ? 

Wave intensity 
ratios preserved? 

Vjli 1 = 

1 

400 

yes 

no 

ambient 

2 

400 

yes 

Yes 

ambient 

3 

400 

no 

Yes 

ambient 

4 

2000 

yes 

Yes 

(ty0x + ambient) /2 

5* 

2000 

yes 

Yes 

ambient 

6 

400 

yes 

Yes 

Vjl 


reference case 





Energy (MeV) r ( solar radii ) 




FIGURE 1 . Evolution of: (a) j'e vs E proton intensity spectrum and (b) I R+ vs k/B Alfv6n wave 
spectrum at the first US cell; (c) j'e vs r profile and (d) mean-free-path X vs r profile of 1.3 MeV protons. 
Panels (e) and (f) are as in (c) and (d) but for 67.5 MeV protons. 


The evolution of I R+ and 1 L ~ DS is illustrated in Figs. 2a and 2b, respectively. These 
wave intensities have been enhanced by the shock transmission process and subject to 
wave saturation while maintaining the ratio of I L ~ (k) / I R+ (k) obtained from (eq. (7)). 
Both DS wave spectra are shifted to higher k out of resonance with the currently highest 
energy protons except at small jl. Furthermore, the I R+ spectrum is shifted more than the 
I L ~ spectrum, so that I L ~ dominate at low k. Consequently the wave compression ratio 
experienced by the currently highest energy protons is weaker than that felt by lower 
energy protons. 

For the reference case, we have employed 200 DS cells and 2000 US cells, ambient 
wave intensities at the moving outer boundary, and preserve the ratios of inward to 
outward wave intensities in applying ad hoc wave saturation. The consequences of 
varying these parameters or procedures on the SEP j’e versus E spectrum at t = 297 s 
are shown in Fig. 2c. The six model runs for comparison are summarized in Table 1 and 
spectrum 5 is for the reference case above. 

Reducing the number of US cells from 2000 to 400 alone produces spectrum 2, which 
shows significant acceleration to only 30 MeV. This can be understood as follows. The 
SEPs amplify the R+, L+ ambient waves as they travel from the outer boundary to the 
shock. The proton intensity scale length is much longer at > 30 MeV than at < 5 MeV 
(Figs, lc and le). Hence, for < 5 MeV protons the resonant wave growth occurs mostly 
near the shock, but for > 30 MeV protons, significant wave growth occurs beyond 
400 cells. Therefore reduction of the US region has little effect on SEP acceleration 










FIGURE 2. (a) DS l L ~ , (b) DS I R+ , (c) je shock spectra at t = 297 s in various simulations (Table 1) 


at < 5 MeV but significantly curtails it at > 30 MeV. Interestingly, accounting crudely 
for wave growth beyond cell 400 using the boundary condition 'F|y= 40 i = 'P/= 4 oo yields 
spectrum 6, which slightly exceeds the reference spectrum 5 at high energies. 

Setting Dfip = Dp M = Dpp — 0 thereby turning off energy loss due to scattering 
by DS inward waves results in spectrum 3. Comparison of spectra 2 and 3 show that 
scattering by DS inward waves reduces SEP acceleration by ~ 30 MeV. Spectrum 1 
results from simply imposing wave saturation without maintaining the ratio of inward to 
outward DS wave intensities predicted by shock transmission. The consequent dominant 
DS R- and L- waves return particles to the shock with energy loss thereby slowing 
acceleration. Finally, spectrum 4 is produced with 2000 US cells and the modified outer 
boundary condition 'P|;=200i = ( v P|y=2000 + ^ambient)/^- Spectra 4 and 5 are almost 
indistinguishable. 


DISCUSSION AND SUMMARY 

Shock transmission has subtly different implications on steady-state and time-dependent 
shock acceleration. Since c w > Cf, steady-state test-particle theory immediately gives 
harder than classical SEP spectral index [6]. In our time-dependent shock acceleration 
model, US SEPs at E < E max amplify outward Alfven waves at k > k m i„. The amplified 
US waves are then shock-transmitted DS with k mm increased more for the outward than 
the inward daughter mode, resulting in (a) loss of resonance with the highest energy 
SEPs except at small jJ and (b) a weaker c w compression at high energy with the DS 
inward mode. This situation is compounded by the SEPs damping (amplifying) DS 
outward (inward) waves, since the DS inward waves return SEPs to the shock with 
energy loss compared to outward waves. Acceleration and wave growth still proceed 
quickly, however, because high-energy SEPs of smaller [i are scattered by these waves. 
These processes are studied via the comparison runs reported above. 

It is important to have q -dependence in the wave-particle resonance condition [12]. 
Quasilinear theory [13] gives B/k ~ P(p - V A /v), so that waves excited by, e.g., 0.5 
MeV protons at j u ~ 0.8 can scatter 2 MeV protons at /r ~ 0.4. This coupling of 
low and high energy protons via the same resonant waves is essential for bootstrap 
shock acceleration and wave growth. If the scattering rate uses the popular “sharpened” 




resonance Bjk — P and the wave growth rate uses the correct /i -dependence [e.g., 14], 
wave growth does not produce commensurate enhancement in the scattering rate for the 
most energetic SEPs and acceleration is slow. Going one step further in the sharpened 
resonance approximation and evaluating the wave growth rate from particles at a single 
energy, much faster acceleration is obtained [15]. However, the pair of wave growth and 
particle scattering rates, based on one-on-one relation between wavenumber and particle 
energy regardless of pitch angle, represents a fairly drastic approximation to the original 
quasilinear expressions. The sharpened resonance prohibits coupling between SEPs of 
different energies via wave interaction, so that in a time-dependent model, the few high 
energy SEPs must amplify their own resonant ambient waves before they experience 
enhanced scattering, slowing down acceleration. 

Self-amplification of upstream Alfven wave turbulence is a powerful process that 
rapidly bootstraps shock acceleration of SEPs. There is likely to be more scattering than 
imposed by the ad hoc A > 3 r g limit in our model [16]. Hence a > 2000 km/s coronal 
shocks is likely to accelerate protons to 100 MeV in < 5 minutes. 

The present model does not consider compressive wave modes, wave-wave inter- 
action, nonlinear cascading, wave dissipation, shock obliquity, or the evolution of the 
shock structure by the SEPs and waves. Extending the model to include any of these 
poses a significant challenge. In focusing on the self-consistent interaction between 
SEPs and Alfven waves, the present work has identified the importance of the interplay 
between SEP-Alfven-wave interaction, shock transmission, and shock acceleration, in 
particular, the role of the downstream waves. 
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